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Generating coefficients for a prediction filter in an encoder 



The invention relates to a transmitter for transmitting a digital information 



signal via a transmission medium, including: 

- input means for receiving the digital information signal, 

- adaptive prediction filter means adapted to derive a prediction signal from the digital 
5 information signal in dependence on an array of prediction filter coefficients, 

- first signal combination means for combining the digital information signal and said 
prediction signal so as to obtain a residual signal, 

- encoding means for encoding said residual signal so as to obtain an encoded signal, 

- coefficient generator means for generating an array of filter coefficients A[i] in response to 
10 the digital information signal, i being an integer for which it holds that 0 < i < p, where p is a 

variable, 

- output means for supplying the encoded signal to an output terminal for transmission via a 
transmission medium. 



15 and generating digital information therefrom, and to a transmission method. 

A transmitting device and a receiving device of the type defined in the opening 
paragraph are known from J. Audio Eng. Soc, Vol. 44, No. 9, pp. 706 - 719, September 

20 1996. The known transmitting device is intended to efficiently reduce the bit rate for the 

transmission of a digital information signal. Prior to encoding the digital information signal, a 
predicted version of the digital information signal is subtracted from the digital information 
signal. The residual signal thus obtained is successively encoded in the encoder and 
transmitted via the transmission medium. The performance of the linear prediction filter is 

25 crucial for the coding gain of the encoder. The performance of the linear prediction filter is 
determined by the prediction filter coefficient. A common method of finding the prediction 
coefficients (a-parameters) is the auto-correlation method. It appeared that a-parameters 
determined with the auto correlation method did not always result in the optimum coding 
gain. 



The invention further relates to a receiver for receiving a transmission signal 
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It is an object of the invention to provide a transmitting device and a receiving 
device using a more efficient method of transmitting and receiving the digital information. 
5 " The transmitting device in accordance with the invention is characterized in 

that the device further comprises smoothing means for smoothing the array of filter 
coefficients A[i] so as to obtain the array of prediction filter coefficients for supply to the 
adaptive prediction filter means. 

The receiver in accordance with the invention is characterized in that the 
10 receiver further comprises smoothing means for smoothing the array of filter coefficients A[i] 
so as to obtain the array of prediction filter coefficients for supply to the adaptive prediction 
filter means. 

The invention is based on the following recognition. As long as transmission 
bandwidth and storage capacity is limited, there will be a need to increase the bit reduction of 

15 digital signals. It appeared that the acquired prediction filter coefficients (a-parameter) 

obtained by means of, for example, the auto-correlation method, did not result in an optimal 
coding gain. During extensive searches on the a-parameter values, it appeared that better 
coding gains could be achieved by slightly modified a-parameters. Coding gain 
improvements of about 4% could be achieved. However, the complexity needed in the 

20 transmitting device for such an extensive search would be enormous, and thus not practical. 

After investigation of the differences between the a-parameters before and 
after the extensive search, it appeared that the a-parameters are roughly the same (see Figure 
8). When zooming in at the first few a-parameters the a-parameters, obtained by means of the 
auto-correlation method show a rippling behavior while the a-parameters resulting from the 

25 extensive search are a smoothed version of it (see Figure 9). 

Instead of executing an extensive search on the a-parameters so as to obtain 
the prediction filter coefficients, the a-parameters acquired by means of the auto-correlation 
method are post-processed by way of a smoothing function. The improvement of the coding 
gain obtained with the smoothing function is comparable with the coding gain obtained with 

30 extensive search. Figure 10 shows the results for a 5-minute DSD fragment, where the result 
was even an improvement of 4.5%. Figure 1 1 shows that the performance is comparable in 
both methods (the extensive search would gain another 0.4% for this fragment). The increase 
of the complexity of the arrangement to implement the smoothing method is negligible as 
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compared with the increase of the complexity of the arrangement when implementing the 
extensive search method. 

These and other aspects of the invention are apparent from and will be 
elucidated with reference to the embodiments described hereinafter. In the drawings, 



value of max |a,-|, 

Figure 4 shows the coding gain as a function of max |a,|, 

Figure 4a shows a detailed view of the coding gain as a function of max |a,| 

shown in Figure 4, 

Figures 5 and 5a show the coding gain and impulse response as a function of 
the prediction order, 

Figures 6 and 6a show the coding gain and impulse response as a function of 
the prediction order after applying the first two steps of the method according to the 
invention, 



coefficients versus the coding gain obtained by using the smoothed version of the first array 
of coefficients, 



obtained by extensive search versus the coding gain obtained by using the smoothed version 
of the first array of coefficients, 



Figure 1 shows an embodiment of an SACD encoder, 
Figure 2 shows the coding gain of successive frames, 

Figure 3 shows the coding gain of successive frames and the corresponding 




Figure 1 1 shows the coding gain obtained by using the array of coefficients 
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Figure 12 shows an embodiment of a transmitter according to the invention, 
Figure 13 shows an embodiment of a receiver according to the invention, and 
Figure 14 shows a flow diagram of an embodiment of a smoothing algorithm. 



For "Super Audio CD" (SACD), the audio signals are lossless coded, using 



framing, linear prediction and entropy coding. The coding gain tj is heavily determined by 
the quality of prediction. It is therefore important to have an algorithm, which optimally 
determines the prediction coefficients. Figure 1 shows a block diagram of a simplified 

10 scheme used in a SACD encoder. The DSD signal, which is a binary representation of an 
audio signal is received at the input terminal 2. The DSD signal is supplied to a coefficient 
generator unit 4 and an adaptive prediction filter 6. The coefficient generator unit 4 is 
arranged to derive coefficient for the adaptive filter such that the prediction filter generates a 
good prediction of the received DSD signal. A common method of finding the prediction 

15 coefficients is the auto-correlation method. The adaptive prediction filter 6 is arranged to 

generate a predicted DSD signal from the received DSD signal. The filter characteristic of the 
prediction filter depends on the supplied coefficients. Signal combination unit 8 subtracts the 
predicted DSD signal from the received DSD signal so as to obtain a residual signal. The 
residual signal is supplied to an encoding unit 10 so a to obtain an encoded signal. The 

20 encoded signal, which is a bit stream signal is supplied to a transmission medium. Depending 
on the implementation of the receiver the coefficients of the adaptive prediction filter have to 
be transmitted. The transmission medium could be a recording medium, such as an optical 
recording medium. The encoding unit comprises preferably a lossless encoder. Such an 
lossless encoder is described in F. Bruekers et al, "Improved lossless coding of 1-bit audio 

25 signals", presented at 103rd Convention of the AES, Sept., 26-29, 1997, preprint 4563(1-6). 
However other coder such as arithmetic coders, huffmann coder are suitable as well. The 
DSD signal is preferably divided into frames. Each frame is applied to the encoder resulting 
in an encoded signal and coefficients which are submitted to the transmission medium. 



suffer from significantly low gains. Such a low coding gain on frames reduces the average 
coding gain and consequently has an adverse implication on the storage capacity. Figure 2 
shows the coding gain of successive frames illustrating the problem. In an excerpt, when the 



Disclosed is a low-complexity method of generating prediction coefficients 
which result in an improved coding gain. 

The problem is that frames, for which high coding gain can be expected, often 
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coding gain is relatively high (average 2.7), it can clearly be seen that there are frames with a 
very low coding gain (close to 1.0). 



problem. However, an important issue is that coding gain is only available once the entire 
5 ~ coding process has taken place and the complexity of interactively/iteratively re-encoding the 
frame using lower order is high. 



coefficients rendering low coding gain, where higher gain was possible, prior to encoding of 
the entire frame. The detection is based on the dynamic range of the direct form prediction- 

10 filter coefficients, often referred to as a-parameters. Consider Figure 3, where the top trace is 
the same as Figure 2 and the lower trace shows the value of the maximum of the absolute 
values of the a-parameters. It is clear from this Figure that the occurrence of low coding gain 
has a correlation with the occurrence of high values of max Based on this detection, a 3- 
stage process is disclosed to eliminate this problem. The result is prediction coefficients 

15 offering higher performance of the encoding process. 



of the absolute values of a-parameters, max \a\, is shown on the horizontal axis and the 
vertical axis shows the corresponding coding gain. The a-parameters are determined by using 
the auto-correlation method in combination with the Schur algorithm. A line is drawn at a 
20 value of 10. Note that for cases, where the maximum of the absolute value of a-parameters is 
higher than 80, there is not a single occurrence of a frame with high coding gain. 



gain and looking at the prediction gain as a function of the prediction order, we see Figure 5. 
Figures 5 and 5a highlight two issues. Figure 5 shows the coding gain as a function of the 

25 prediction order. Figure 5a shows the impulse response for various prediction orders. One 

may conclude from Figure 5 that a higher coding gain was possible, for instance at prediction 
order 80, where the gain was maximum. The second issue is highlighted in Figure 5a. The 
impulse response only changes in scale but not in shape. It is interesting to interpret this in 
the spectral domain, if one realizes that the prediction gain is proportional to the non-flatness 

30 in the spectrum [Flanagan paper]. Thus, all the flattening has already been achieved in lower 
order, and when the order is increased, only the scale changes while the spectral coloring (or 
the flattening in the inverse filtered domain) does not change. This points to some ill 
conditioning. Steps were taken to avoid this problem by inspecting the impulse response of 
the prediction filter, i.e. without having to go through the entire coding process. 



Interactively, selecting a lower prediction order would have prevented this 



A detection mechanism is disclosed, which is capable of signaling LPC 



Figure 4 illustrates a slightly different version of Figure 3. Here, the maximum 



When zooming in to one of the "problem frames", i.e. a frame with low coding 
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There appeared to be a relation between low coding gains and the maximum 



absolute value of the a-parameters returned by the auto-correlation method. For frames with a 
- good coding gain, the value of the maximum absolute a-parameters is always below 10, and 
problem frames show values of 100 or even 600. 



the problem frames, i.e., frames for which the coding gain is low, while a higher gain was 
possible. Accordingly, step 1 includes the computation of the a-parameters based on the auto- 
correlation method. Next, a value T, representing the absolute maximum value of a, is 
computed on the basis of the a-parameters, A = { aj f ci2,..., am h such that, T = max |A|. A 

10 frame is declared as a problem frame if T is larger than a predefined value, for example 10.0. 

Step 2: Some of these problems are eliminated if one adds a very low-level 
random noise signal to the input prior to LPC analysis. Equivalently, the value of R[0], the 
first auto-correlation coefficient, can be raised by a very small amount. If the threshold of 
step 1 is exceeded, the auto-correlation function coefficient R[0] is modified in such a way 

15 R'[0] = R[0] * (1.0 + 4.10" 6 ), and the a-parameters are recalculated in accordance with the 
new value of R[0], 

Although such a minor modification solves the problem in most cases, there 
were frames for which the problem persisted, or even got worse (in terms of coding gain). 

Step 3: The maximum absolute value of the now acquired a-parameters is 
20 again a good detection mechanism to see if problems are to be expected. For these frames, 
where the addition of a little noise does not solve the problem, another measure must be 
taken and that is to take a reduced order for LPC prediction. It is difficult, however, to select 
the optimal prediction order. A threshold which seems to achieve good results is: 
New Order = min ( (current order * 8 / 9) , 80 ) where the current order is the order of 
25 prediction with which the analysis has been done. 

In practice it appeared that with the combination of these three steps (detection 
via max(abs(a-parameters)), noise addition, and order reduction), the coding gain improved 
considerably for problem frames. More importantly, the good frames were not affected 
adversely with this mechanism. 
30 Results: 

Figure 6 shows the improvements achieved after applying the first 2 steps of 
the procedure discussed above. This is in contrast to Figure 5 where the coding gain dropped 
sharply beyond 80 th -order prediction. Now the a-parameters also exhibit a reasonable 



The step 1 in the method in accordance with the invention is the detection of 
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behavior in terms of the amplitude dynamics. As indicated above, even after applying the 2 
steps, some problem frames persist and the application of step 3 eliminates them as well. 

Finally, Figures 7, 7a and 7b show clearly that the bad frames were completely 

eliminated. 

A three-stage approach of eliminating the bad frames entirely has been 
described. This was achieved by conditioning the signal in such a way that prediction 
coefficients are obtained which result in an improved coding gain. 

While only Schur recursion has been discussed, experiments show that 
identical problems exist when using the Cholesky decomposition. Also in this case, the 
disclosed solution helps to completely eliminate the ill-conditioned a-parameters. 

"Smoothing" of the a-parameters or prediction filter coefficients could yield 
(about) the same results as the extensive search, through for a negligible complexity. 

A smoothing algorithm uses the a-parameters for calculating a new set of a- 
parameters. Possible implementations could be done by, for example, FIR or IER-filtering. 
After some experiments, a variant of the IIR-method was selected: 

/* Apar[0] .. Apar[po-l] contain the a-parameters of the filter with prediction order po */ 
for (i = 1; i < po - 1; i++) { Apar[i] = Apar[i-1] + 2*Apar[i] + Apar[i+1] } 

As can be seen from this formula, the algorithm is recursive because the value 
which is adapted is reused in the next calculation. It appeared that it had a better effect than 
the non-recursive approach. 

On average, some 3% coding gain improvement can be reached when 
applying this smoothing to the a-parameters. Figure 10 shows the results for a 5-minute DSD 
fragment, where the result was even an improvement of 4.5%. The complexity of this method 
is negligible. As can be seen from Figure 10 there were no frames where smoothing 
decreased the coding gain. Comparable results were obtained for some other DSD-fragments, 
with the exception that the gain showed a minor decrease for a very few high gain frames 
0l>3). 

A comparison is made between the coding gains resulting from the extensive 
search and those resulting from smoothing the a-parameters that resulted from the auto- 
correlation method. Figure 1 1 shows that the performance is comparable in both methods (the 
extensive search would gain another 0.4% for this fragment). To decrease the complexity of 
the arrangement, it is obvious that the method using smoothing of the a-parameters is 
preferred. 
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Figure 14 shows a block diagram of another embodiment of the smoothing 
algorithm. The smoothing algorithm generates n coefficients for the prediction filter having n 
. coefficients. Ci n [0:n-1] are the n filter coefficients before smoothing. C ou t[0:n-1] are the filter 

coefficients after smoothing. The following algorithm obtains the smoothed coefficients: 
" CouttO] = Q n [0], 
C out [n-l] = Q n [n-l], 

C ou t[i] = 0.25*C in [i+l] + 0.5*C in [i] + 0.25*C out [i-l], 
whereby I is an integer and 1 < i < n-2. 

Figure 14 shows the block diagram representing the last expression. The 
square blocks represent a delay of one coefficient. The triangles represent the multiplication 
of the input signal with the corresponding factor. The ellipse represents the addition of all the 
inputted signals. 

Figure 12 shows an embodiment of an encoder in accordance to the invention. 
The encoder of Figure 12 shows similarity with the encoder of Figure 1. The units having the 
same reference number perform the same function. The encoder shown in Figure 12 differs 
from the encoder shown in Figure 1 in that the linear prediction coefficients are applied to a 
smoothing unit 12 prior to supplying the coefficients to the prediction filter 6. An 
embodiment of the smoothing function performed by the smoothing unit is described above. 
Preferably the smoothing unit performs a low-pass filtering on the coefficients. Such filters 
could be implemented by finite impulse response filters (FIR), infinite response filter (IIR) or 
a combination of FIR and IIR. Figure 14 shows a combination of an IIR and a FDR filter. 

Figure 13 shows an embodiment of a receiver in accordance to the invention. 
A transmission signal comprising an encoded signal, is received at an input terminal 18. The 
received signal is supplied to a decoder unit 20. The decoder unit 20 decodes the encoded 
signal so as to obtain a residual signal. The decoder unit preferably comprises a lossless 
decoder. Such an lossless decoder is described in F.Bruekers et al, "Improved lossless coding 
of 1-bit audio signals", presented at 103rd Convention of the AES, Sept., 26-29, 1997, 
preprint 4563(1-6). However other coder such as arithmetic decoders, huffmann decoder are 
suitable as well. The residual signal is supplied to a signal combination unit 22. The signal 
combination unit is arranged to add a predicted digital information signal to the residual 
signal so as to obtain the digital information signal. The digital information signal is supplied 
to an output terminal of the receiver 30. The digital information signal is further supplied to a 
coefficient generator unit 26 and an adaptive prediction filter 24. The coefficient generator 
unit is arranged to derive filter coefficients for the adaptive prediction filter from the digital 
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information signal. A common method of finding the prediction coefficients is the auto- 
correlation method. However others method are suitable as well. The prediction filter 
- coefficients are supplied to a smoothing unit 28 adapted to perform a smoothing function on 
the prediction coefficients so as to obtain smoothed prediction filter coefficients. Preferably 
the smoothing unit performs a low-pass filtering. The smoothed prediction filter coefficients 
are supplied to the adaptive prediction filter 24. The adaptive prediction filter is adapted to 
generate the predicted digital information signal. It should be noted that for lossless 
transmission of the digital information signal the adaptive prediction filter, the coefficient 
generator unit and the smoothing unit in the transmitter and receiver perform exactly the 
same functions. It should further be noted that, the coefficient generator unit 26 and the 
smoothing unit 28 are optionally if the transmission signal comprises an encoded signal and 
the adaptive prediction filter coefficients. In said case the adaptive prediction filter 
coefficients need not to be derived by the receiver from the digital information signal. 



alternative that where entire coding process has to be performed before the coding gain is 
known. 



thereof, it is to be understood that these are non-limitative examples. Thus, various 
modifications are conceivable to those skilled in the art, without departing from the scope of 
the invention, as defined by the claims. 



article "a" or "an" preceding an element does not exclude the presence of a plurality of such 
elements. In the claims, any reference signs placed between parenthesis shall not be 
construed as limiting the scope of the claims. The invention may be implemented by means 
of hardware as well as software. Several "means" may be represented by the same item of 
hardware. Furthermore, the invention resides in each and every novel feature or combination 
of features. 



The disclosed method is also lower in complexity when compared with the 



Though the invention is described with reference to preferred embodiments 



The use of the verb "to comprise" and its conjugations does not exclude the 
presence of elements or steps other than those stated in a claim. Furthermore, the use of the 



